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Abstract: We consider the problem of implementing mutually unbiased 

bases (MUB) for a polarization qubit with only one wave plate, the mini¬ 
mum number of wave plates. We show that one wave plate is sufficient to 
realize two MUB as long as its phase shift (modulo 360°) ranges between 
45° and 315°. It can realize three MUB (a complete set of MUB for a 
qubit) if the phase shift of the wave plate is within [111.5°, 141.7°] or its 
symmetric range with respect to 180°. The systematic error of the realized 
MUB using a third-wave plate (TWP) with 120° phase is calculated to be 
a half of that using the combination of a quarter-wave plate (QWP) and 
a half-wave plate (HWP). As experimental applications, TWPs are used 
in single-qubit and two-qubit quantum state tomography experiments and 
the results show a systematic error reduction by 50%. This technique not 
only saves one wave plate but also reduces the systematic error, which can 
be applied to quantum state tomography and other applications involving 
MUB. The proposed TWP may become a useful instrument in optical 
experiments, replacing multiple elements like QWP and HWP. 
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1. Introduction 

Measurement bases with a special geometry structure in Hilbert space, such as 2-designs m, 
weighted 2-designs a, tight frames II and Platonic solids a, have attracted wide attention 
in the community of quantum information science in recent years. A typical example of 2- 
designs is mutually unbiased bases (MUB) @ 13 . The special geometry for MUB is the equal 
Hilbert-Schmidt overlap between any two projective operators corresponding to two MUB. 



Hence, measurement results of one operator provide no information about the measurement 
results of its MUB at all. Due to this geometry, MUB have found applications from quantum 
key distribution lT8ffl~Qjl . entropy uncertainty relations mm to quantum state estimation 10 

mna. 

A complete set of MUB in a multi-qubit system involves nonlocal measurements that are 
difficult to be realized in the laboratory. In many scenarios such as quantum state tomography 
(USED, the tensor product of a complete set of single-qubit MUB is a preferred choice as 
the measurement bases for multi-qubit systems mini. Therefore, it is important to implement 
MUB for a single qubit in experiments ll22U24ll . 

In the context of polarization optics, a specific type of polarization transformation poses 
constraints on the minimum number of optical elements that must be used. Three wave plates 
are the minimum number to realize any SU(2) polarization transformations ll25ll and the visual 
tool kit can be found in [261. As for transformations between a given pair of nonorthogonal 
polarization states, two QWPs are enough (see a sketched proof in [27] and a recent constructive 
demonstration by Zela in |28| ). When the initial state is linearly polarized, a quarter-wave plate 
(QWP) and a half-wave plate (HWP) can realize any polarization state ED . The combination of 
QWP and HWP is used in most of current optical experiments f29l where the linear polarization 
is generated by a polarizing beam splitter (PBS). In the setting composed of a QWP, a HWP and 
a PBS, the PBS acts as <7- and by adjusting the optic axis angles of the QWP and HWP (a two- 
parameter setting), the Bloch vector of <7- can be unitarily rotated to any point on the unit Bloch 
sphere ED- However, when the problem is restricted to realize MUB rather than arbitrary 
bases, is one wave plate sufficient to transform the initial <7- to a set of MUB, especially a 
complete set of three MUB? If the answer is positive, not only one wave plate is saved but also 
the systematic error H29I130I can be expected to decrease as the parameter uncertainties of the 
devices (i.e., wave plates in our case) are the sources of systematic error. Intuitively, the less 
measurement devices, the smaller the systematic error. So far, people know how to construct 
two MUB with one HWP or one QWP. For example, one HWP itself can be used to realize el¬ 
and o x by setting its rotation angles as 0° and 22.5°; one QWP can be used to realize O- and 
O y by setting its rotation angles as 0° and 45°. One wave plate is also used to perform Fourier 
transform tomography |3J) and the optimal phase of the wave plate is numerically calculated 
in (32l . As to a complete set of three MUB, this problem is barely considered. In this paper, 
we show that it is indeed feasible to employ only one wave plate to realize a complete set of 
single-qubit MUB with reduced systematic error. 

Here is the organization for the rest of the paper. In section [2j the conditions to realize two 
MUB and three MUB using only one wave plate are considered, respectively. Section [3] cal¬ 
culates the systematic error in the realization of MUB with only one wave plate. In section 
[4] third-wave plates are used to perform MUB measurements in both single qubit and two- 
qubit tomography experiments and the results demonstrate an error reduction by 50%. Section 
[^concludes the paper. 

2. Theory 

2.1. Transforming polarization states with QWP-HWP setting 

In optical experiments, arbitrary projective measurements on polarization qubits are often im¬ 
plemented by a QWP-HWP setting; see Fig. |TJa). Define | H) = (l,0) r , |V) = (0, l) r and their 
corresponding eigenvalues are ±1, where T denotes transpose. Then the PBS acts as (7- or 
(0,0,1 ) 7 in the Bloch representation. This setup transforms | H) to 
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QWP HWP PBS WP PBS 



(a) (b) 


Fig. 1. MUB measurement with two different settings: (a) MUB measurement with QWP- 
HWP setting consisting of a QWP, a HWP and a PBS. (b) MUB measurement with a wave 
plate (WP) and a PBS. 


where qji are, respectively, the optic axis angles of QWP and HWP deviated from the horizon¬ 
tal direction, and U is a unitary transformation operator on polarization by a phase plate with 
phase 5 and rotation angle 0, 


U(0,8) 


cos 2 9 + e iS sin 2 0 |(1 - e iS ) sin 2 6 \ 

|(1 — e' 5 )sin20 sin 2 9 + e ,s cos 2 0 )' 


( 2 ) 


Or equivalently this setup transforms the original operator <7- to the Pauli operator r(q.h) ■ a 
with 

r(q,h) = (sm2qcos(4h — 2q),sm(4h — 2q),cos2qcos(4h — 2q)) T (3) 

and a = (oy, ay, a,), where ay, ay and o- are three Pauli operators. 

Two variables of rotation angles q and h correspond to a plane and can rotate the initial Bloch 
vector (0,0, l) r to any direction or position on the unit Bloch sphere. Thus, this QWP-HWP set¬ 
ting is used to realize any one qubit projective measurement basis. For example, the most popu¬ 
lar set of MUB is composed of Bloch vectors (l,0,0) r , (0, l,0) r and (0,0, l) r (corresponding 
to three Pauli operators) which are realized by choosing q = 45°,0°,0° and h = 22.5°,22.5°,0° 
correspondingly. 

In order to realize MUB for a qubit, only two or three specific Bloch vectors need to be 
implemented and it is not necessary to cover the Bloch sphere with two rotation angles of two 
wave plates. We will show below that one wave plate with appropriate phases is sufficient to 
realize single-qubit MUB. 


2.2. Realizing two mutually unbiased bases with one wave plate 

Two sets of orthogonal bases {|y / /)}fLi and {| )}y=i i n a ^-dimensional Hilbert space are 

mutually unbiased if |(V / , 1 |V / /)| 2 = ^ for any i and j. There are at most d + 1 bases with ev¬ 
ery two of them mutually unbiased, which are called a complete set of MUB. So far, how to 
construct a complete set of MUB is known only in systems with dimensions which are powers 
of primes 0E). For general cases, the existence of such a set is still an open problem even in 
the simple case of d = 6 @1121. For a qubit with dimension d = 2, the bases {|t/^/)} 2 = i are 
the eigenvectors of a unit Pauli operator r 7 ■ <7 with | r J \ j = 1. In the Bloch representation, the 
yth basis is directly related to its Bloch vector r J as i)//)( \j//\ = I±r ^' a and ± corresponds to 
i— 1,2. Thus, the requirement for two bases to be mutually unbiased is that their Bloch vectors 
are normal to each other 

r 1 • r 2 = 0. (4) 

In Fig. Hlb), a wave plate (WP) with a phase difference 8 between ordinary (o) and ex¬ 
traordinary (e) components is placed in front of a PBS. Its unitary transformation operator on 






















Fig. 2. Frame potential *5(5) with respect to 8 . The frame potential is symmetric about 
180°. The inset is the amplified interval of 8 where the frame potential vanishes to zero. 
Only if 8 is in the range [111.5°, 141.7°] or its symmetric range about 180° can a complete 
set of MUB be realized by one wave plate. 


polarization is expressed in Eq. (0 where 0 is the deviation angle of the optic axis from the 
horizontal direction. This wave plate with rotation angle 9 transforms the initial Bloch vector 
ro = ( 0 , 0 , l) r into 

r(9) = (sin 2 fsin40, -sin<5sin20, 1-2sin 2 §sin 2 20) r . (5) 

From Eq. ©, this setting with one wave plate can realize two MUB if there exist 9\ and 02 
such that their Bloch vectors are normal to each other, i.e., r(0i) • r( 02 ) = 0. The existence 
of 0] and 02 obviously depends on 5. One trivial example is 5 = 0° where the wave plate 
does not affect the initial polarization and r(0i ) • r( 02 ) = 1 for any 9 \, 02 - Another example is 
8 = 180°. The choice of 0j = 0° and 0 2 = 22.5° makes r^) = (0,0, l) r , r(0 2 ) = (1,0,0) T , 
and thus r(9\) ■ r{9i) = 0. The maximum of r(9 1 ) ■ r( 02 ) is 1, which can be achieved at 9\ = 02 - 
Because r(9 \) ■ r( 02 ) is a continuous function and its maximum is positive, there exist 0] and 
02 to make r(0i ) ■ /"(Oo) = 0 if the minimum of r(0i ) • r( 0?) is negative or zero. This problem is 
then converted into finding the minimum of r(0i ) • r(0o). From Eq. 0. the first two elements 
of r are odd functions of 0 and the third term is an even function of 0. That is, n(— 0) = 
—n( 0 ),r 2 (— 0) = —^ 2 (0) and r 3 (—0) = 7 - 3 ( 0 ). From Eq. 0. the minimum and maximum 
of 7 - 3 ( 0 ) is cos 8 and 1. If cos 5 is non-positive, then there exists 9\ = ^ arc sin (^ esc |) that 
makes r 3 (0i ) = 0 and we choose 02 = — 9\ to make r(0i) = — r( 02 ). We obtain the minimum 
r(0i) -r( 02 ) = — 1. Thus, in this case, two MUB can be realized such as 0i = 7 arcsin(-^csc 4 ) 
and 02 = 0°. If cos 8 is positive, the minimum of r (0i) • r ( 02 ) is cos 2 5 obtained at 0i = — 02 = 
45°. Thus, the requirements that the minimum of r(0i) • r( 02 ) is non-positive in these two cases 
together give 45° < 8 < 315°. As long as the phase 5 (modula 360°) of the wave plate is within 
[45°,315°], it can be used to realize at least two MUB. As there are two parameters and one 

equation, the solutions of Eq. 0 form a line. One solution is 0] = ^arcsinf ^ cscl)) and 

02 = —0i, which make 7-3 = The Bloch vectors of the two MUB are axially symmetric 
about the initial state (0,0, l) r , and the angles between MUB and (0,0,1 ) 7 are 45°. 










180 



Fig. 3. The numerical solutions of rotation angles in Eq. © with respect to the phase 8. 
There are two classes of solutions with S between 111.5° and 126.3°, represented by solid 
and dotted lines, respectively, and only one class of solutions in (126.3°, 141.7°]. Each 
class of solutions contains four complete sets of MUB, denoted by four different colors. 


2.3. Realizing a complete set of MUB with one wave plate 

Three rotation angles 0,-, i = 1,2,3, are chosen to realize three MUB, which consist of a com¬ 
plete set of MUB for a qubit, 

r(6i)-r(dj)= 0 , ( 6 ) 

where i. j = 1,2,3 and i / j. For a complete set of MUB, it is much more complicated and 
difficult to obtain theoretical solutions. Hence, the problem is investigated numerically in this 
paper. Borrowed from frame theory [f2l|34| , the frame potential of a wave plate is defined as 

0(5) = min [r( 0 [) ■ r( 0 2 )] 2 + [r( 0 2 ) ■ r( 0 3 )] 2 + [r( 0 !) ■ r( 0 3 )] 2 . (7) 

0 ]. 02 5 03 

The solution of Eq. 0 exists only if the phase 5 of the wave plate makes the frame poten¬ 
tial Eq. £7} vanishing. The frame potential is numerically computed by a MATLAB solver 
Isqnonlin, which is intended to solve nonlinear least-squares problems. One hundred different 
initial points of 6\, 0 2 and 0 3 are taken for us to avoid local minimum and find all the pos¬ 
sible solutions. The numerical result of the frame potential for different phases is shown in 
Fig. [2] From Eq. (0, $(5) has a period of 360°. Thus, we only consider 0° < 8 < 360°. Since 
r(0,,5) ■ r(9j, 8) = r(0,-, 360° — 8) ■ r{9j, 360° — 8), we have 0(5) = 0(360° — 8), meaning 
that the frame potential 0(5) is symmetric about 5 = 180°. 

With 5 between 111.5° and 141.7°, the frame potential is zero and the solutions of Eq. 0 
are shown in Fig. [3] From Eq. 0, r(0) = r(9 + 180°), we only consider rotation angles within 
[0°, 180°). Eq. © also shows r(0 ; ) • r(0,) = r{90° ± 0,-) • r(90°± 9j) = r(180° - 0,) ■ r(180°- 
9j). That is to say, if 0|. 0 2 . 0 3 (red in Fig. 0 is the solution, the modules of 90° ± 0,- and 
180° — 0, (i = 1,2,3) by 180° (represented as blue, green and black in Fig0correspondingly) 
are also the solutions. Thus, all the four sets of MUB are considered as one class of solutions. 

As shown in Fig© solutions from different sets of MUB in the same class intersect around 
5 — 120°, 126.3° and 141.7°. Using the symmetries represented by the colors, the four variables 























(i.e. 8, 0i, 02 and @ 3 ) reduce to two and we can theoretically calculate the solutions at these in¬ 
tersections. The phases at these intersections are also rigourously found to be 8 = 120°, 126.32° 
and 141.76° (see Appendix IA. ll and lA.2l >. Here we compare our results with those in (32). The 
optimal phase was numerically calculated in (32) to be ln/10 (i.e. 126°). The figure of merit 
times the total number of counts in (32) at this optimal phase equals 10.03, which is very close to 
10, the bound achieved by MUB. This phase falls within our calculated range [111.5°, 141.7°]. 
The reason why they only found one phase rather than an available interval of phases and the 
optimal performance at this phase was slightly worse than the bound is their restriction of six 
equally spaced rotation angles. 

In the special case of 8 t = 120°, called third-wave plate (TWP), the Bloch vector is 

( 3 V3 3 1 \ 

- sin4f,—— sin2f,-cos4f + - j , (8) 

where t is the rotation angle of the optic axis of TWP deviated from horizontal direction. 
From Appendix lATl solutions of Eq. © are r(0°) = r(90°) = r(180°) = (0,0, l) r ,r(f 0 ) = 
—r(180° -t 0 ) = ^(l,-l,0)V(f 0 + 90°) = -r(90° - to) = ^(1,1,0) r , where t 0 = 

^ arccos(—ss 27.37°. 


3. Systematic error in the realization of MUB in one wave plate setting 

Imperfect measurement devices are the main sources of the systematic error in the realization of 
MUB. Here we consider the systematic error due to the parameter uncertainties of wave plates. 
The realized bases are denoted by their Bloch vectors as r(5, 0), where 8 is the real phase of 
the wave plate in the one wave plate setting in Fig.QJb). The systematic error in the realization 
of r is 

(Ar) 2 =£ll^H 2 ( A ^) 2 > (9) 

where £ = 5,0. From Eq. (©, 


I f 5 ll = skr20 = A 


(ri) 2 


■ 2 r ’ 

sin o 


/S /S 

||r 0 || 2 = 16sin 2 — —4sin 2 5sin 2 20 = 16sin 2 — —4(r 2 ) 2 . 


( 10 ) 


For a complete set of MUB (rL j = 1,2,3), the systematic error sums up to 

e 2 =t (Ar') 2 = £ I K H 2 (A<I) 2 = E(£ ? ) 2 (A^) 2 , 


2=1 

(«%) 2 = LKl | 2 . 

7=1 


7=1 « 




( 11 ) 


As r 1 is orthogonal to each other, from Eq. ( 1 1 01 ). 


te ) 2 = E-^ = -U, 


(ee) 2 = E 16 sin 2 --4(^) 2 =48sin 2 --4. 

7=i 1 L 





Thus, the systematic error in the one wave plate setting is 


e 2 = -^-(A8) 2 + (48sm 2 ^-4)(A9) 2 . (12) 

sin“ 0 2 

As £ 2 in Eq. (fl2l> is an increasing function of 8 in the interval [111.5°, 141.7°], the minimum 
and maximum systematic error in the realization of three MUB with one wave plate setting is 
1.16(A<5) 2 + 28.8O(A0) 2 and 2.6O(A<5) 2 + 38.83(A0) 2 , achieved at 8 = 111.5° and 8 = 141.7°. 
For a third-wave plate with 8t = 120°, 

£ 2 = 1.33(A<5,) 2 + 32(Af) 2 . (13) 

Under the assumption that (A<5/,) 2 = (At),.,) 2 = (A 8 t ) 2 = (At)) 2 and (Ah) 2 = (A q) 2 = (At) 2 = 
(AO) 2 , the systematic error in the realization of MUB is 1.33(At)) 2 + 32(A0) 2 in the TWP 
setting and averaged as (2.5A5) 2 +68(A0) 2 in the QWP-HWP setting in Eq. d26l >. Thus, the 
TWP setting outperforms the QWP-HWP setting by about a factor of two. 

Measurements based on single-qubit MUB are preferable choices in quantum state tomogra¬ 
phy. In qubit state estimation, a complete set of single-qubit MUB is used to extract information 
of the qubit optimally. In multi-qubit quantum state tomography, the product measurements of 
single-qubit MUB on each photon are used to reduce estimation error due to statistical fluctua¬ 
tion. When the copies of states p are infinite, the estimated state p based on the measurement 
data should be the same as the real state p. However, since single-qubit MUB are imperfectly 
realized, the estimated state p no longer converges to the real state p and tr(p — p) 2 is defined 
as the systematic error in state estimation. Generally, tr(p — p) 2 depends on p. Averaged over 
unitarily equivalent states, (tr(p — p) 2 ) in both single-qubit and multi-qubit state estimation is 
proportionate to the systematic error of the realized bases f35) . As the systematic error in the 
realization of multi-qubit product bases is the sum over the systematic error of single MUB for 
each qubit, the systematic error in the realized multi-qubit product bases with one wave plate 
for each qubit is still a half of that with the QWP-HWP combination. This systematic error re¬ 
duction effect in quantum state tomography is experimentally verified in both single-qubit and 
two-qubit tomography experiments in the next section. 

4. Quantum state tomography experiments with third-wave plates 

4.1. Qubit tomography experiments 

The experimental setup, shown in Fig. [4] includes two parts: state preparation and MUB 
measurement. A 40 mW, V-polarized beam at 404 nm from a semiconductor laser pumps a 
type I phase-matched /3-barium borate (BBO) crystal. After the spontaneous parametric down- 
conversion (SPDC) process, a pair of 808 nm H-polarized photons are created. One photon 
passes through a 3 nm interference filter and is detected by a single photon detector to her¬ 
ald the presence of its twin photon. The quantum state of the heralded photon is prepared by 
HWPO, HWP1, QWP1 and a 770A quartz crystal which is much larger than the coherence 
length of about 270A with A = 808 nm. HWPO with rotation angle h {) and the quartz crystal 
with optic axis aligned horizontally together prepare the quantum state p wih Bloch vector 
s = cos4/i()(0.0. l) r ; HWP1 and QWP1 can transform s to arbitrary direction. This part is ca¬ 
pable of preparing arbitrary qubit state. In the MUB measurement part, a complete set of MUB 
is performed with two methods: a TWP with rotation angles set as 0°, 27.37° and 117.37° 
in Fig. Ufa); the conventional QWP-HWP setting with rotation angles set as (45°,22.5°), 
(0°,22.5°) and (0°,0°) in Fig.Hb). 

As the systematic error in quantum state tomography experiments depends on the state to 
be measured, we prepared three states p = p\(j>)(<j>\ + (1 — p)I/2 at p = 0.92 (by setting 




Fig. 4. Experimental setup for qubit tomography. The apparatus consists of two parts: state 
preparation (green) and MUB measurement (pink). The MUB measurement part consists 
of a polarizing beam splitter and a wave-plate combination which has two choices: (a) TWP 
and (b) QWP-HWP combination. 



p = p|D)(D| + (l-p)7/2 





Fig. 5. Systematic error in TWP and QWP-HWP setting for qubit tomography. The depen¬ 
dence of the systematic error in the TWP setting (Fig. |4ja)) and the QWP-HWP setting 
(Fig.[4jb)) is experimentally measured with respect to the angle errors of TWP (red), QWP 
(black) and HWP (blue) for three states in Figs.[5ja)-{5]c) at p = 0.92. The total systematic 
error (green) in the QWP-HWP setting is the sum of that due to QWP and HWP. The ex¬ 
perimental results denoted as dots coincide with the theoretical calculations (solid lines). 
Fig.[5jd) plots the total systematic error for all these three states in the two settings and 
shows that the TWP setting beats the QWP-HWP setting by a factor of two. Error bars are 
the standard deviation of 100 trials in Monte Carlo simulation with binomial distribution 
of counting statistics. 





























































Fig. 6. Experimental setup for two-qubit tomography. A singlet state is prepared via SPDC 
process with a fidelity of 98%. In quantum state tomography, single-qubit MUB measure¬ 
ments are implemented on both photons with TWP or the combination of QWP and HWP 
in Figs.[4](a) and[4](b). Coincidence events are recorded by a coincidence circuit. 


HWPO’s rotation angle ho = 5.8°) and 10) = //), | D) and | R) with a fidelity of more than 0.998. 
The theoretical systematic error tr(p — p ) 2 in the estimation of these three states is derived in 
Appendix[C]as [8(A h) 2 +4(A < 7 ) + 0.25(A5/,) 2 ]p 2 , [2(A q) + 0.25(A5/,) 2 + 0.5 (A8 ci ) 2 ]p 2 and 
[16(A h) 2 +4(A q)]p 2 in theQWP-HWP setting and [8(Af) 2 + \(A8,) 2 ]p 2 , [^(Af) + |(A ^) 2 ]/? 2 
and [|(Af) 2 + g(A 8,)\p 2 in the TWP setting. Under the assumption that (A5/,) 2 = (A5^) 2 = 
(At),) 2 = (At)) 2 and (Ah) 2 = (At/) 2 = (Ar) 2 = (A0) 2 , the total systematic error of these three 
states is [(At)) 2 + 34(A0) 2 ]/? 2 in the QWP-HWP setting and [| (At)) 2 + 16(A 0) 2 ]p 2 in one wave 
plate setting. That is, the systematic error in realizing the corresponding MUB using the TWP 
setting is around a half of that using the QWP-HWP setting. 

As phase errors are determined by the manufacture and wavelength, without a variable wave¬ 
length we can only experimentally verify the relationship between the systematic error of the 
estimated state and the angle errors of QWP, HWP and TWP for all these three states. In the 
experiment, MUB measurements are performed on 3 x 10 6 photons with two different settings 
in Figs. |4t a) and[4](b). We first measure the state with the well-calibrated setting and assume 
the estimated state p as the real state. Then we intentionally mis-calibrate the optic axes of 
the wave plates with an angular error, and obtain an estimation p. Thus, the systematic error 
due to this angular error is calculated as tr(p — p) 2 . The estimated states by the well-calibrated 
TWP setting have a fidelity of over 99.9% with those by the well-calibrated QWP-HWP setting 
for all the three states above, validating each other. In terms of systematic error, experimental 
results (dots) and the theoretical results (solid lines) are shown in Fig. [5] and they match very 
well. The performance of these two settings depend on p and neither always outperforms the 
other. For example, for states in Figs. [5ja) and 0c), TWP beats the QWP-HWP combination 
while reversely for the state in Fig. [3b). However, there are more states where TWP performs 
better. The total systematic error in the estimation of these three states in the TWP setting adds 
up to be about two times smaller than that in the QWP-HWP setting as shown in Fig.0d). 

4.2. Two-qubit tomography experiments 

In Fig. 0 a 100 mW, H-polarized beam at 404 nm from a continuous laser pumps a pair of 
type I phase-matched -barium borate (BBO) crystals whose optic axes are normal to each 
other. After the spontaneous parametric down-conversion (SPDC) process, a pair of 808 nm 
photons are created. When the optic axis of half-wave plate (HWP0) at 404 nm is deviated 
22.5° from horizontal direction, the twin SPDC photons are maximally entangled. HWP1 and 
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Fig. 7. Systematic error in TWP and QWP-HWP setting for two-qubit tomography. In (a), 
systematic errors due to angle errors of QWP (black), HWP (blue), both of them (green) and 
TWP (red) are numerically simulated and denoted as solid lines. The experimental results 
due to angle errors of wave plates are plotted as dots for photon 1 and circles for photon 2. 
The total systematic error for the two photons is plotted in (b). From (b), The TWP setting 
beats the QWP-HWP setting by about a factor of two. Error bars are the standard deviation 
of 100 Monte Carlo simulations with multinomial distribution of the counting statistics. 


HWP2 rotate H and V to the fast and slow axes of the single mode fibers. At the output ports of 
the fibers, HWP3 and HWP4 rotate the polarization direction back to horizontal and vertical. 
QWP1 is tilted to compensate the phase of the entangled states to a singlet state. In the MUB 
measurement part, a complete set of MUB on either photon is performed with two methods: 
TWP setting in Fig. [4} a) and QWP-HWP setting in Fig. |4jb). 

In the two-qubit case, the singlet state p = 'P ) (*P is chosen for three reasons: firstly, 
systematic error in the estimation of product states is a direct sum of that of single-qubit states; 
secondly, entangled states reveal the peculiar features of quantum systems and are valuable 
quantum resources; the last reason is that the systematic error in the estimation of Werner states 
is proportionate to the sum of the systematic error in the realization of single-qubit MUB for 
either photon (35], which is similar to the systematic error averaged over unitarily equivalent 
states. The systematic error for Werner states l35l is 

tr(p-p) 2 = ip 2 ( £ 2 + e 2 2 ), (14) 

where p = p | X P~) ( V F~| + (1 — p)l/ 4, I 1 ? - ) = -±=(\HV) — | VH)) and ef is the systematic error 
of the realized single-qubit MUB on photon i, i = 1,2. 

Similar to the qubit tomography experiment, we only experimentally measure the depen¬ 
dence of the systematic error on the angle errors of wave plates in the branches of both photons. 
Then product measurements of MUB are performed on 9 x 10 5 pairs of prepared singlet states. 
The estimated state by the well-calibrated TWP setting has a fidelity of over 99.8% with that 
by the well-calibrated QWP-HWP setting, agreeing well with each other. From Eqs. O, o 
and (125b - angle errors of wave plates for either photon theoretically contribute 8(Af) 2 in the 
TWP setting and 12(A/;) 2 + 5(A<r/) 2 in the QWP-HWP setting. However, numerical results in 
Fig.[7](a) show that angle errors only cause 1 0.5 (Ah) 2 + 4.6(Aq) 2 in the QWP-HWP setting and 
6.9(A t) 2 in the TWP setting. This gap is due to the positive semi-definite conditions of density 
matrices, which arises when states are singular. Both the experimental results and numerical 
results in Fig.[7](b) show that the systematic error in the TWP setting is only about a half of 












that in the QWP-HWP setting. The experimental results are slightly smaller than the numerical 
results because the prepared state is not exactly the expected singlet state, which only has a 
fidelity of 98%. 

5. Conclusion 

We have found that one wave plate is sufficient to realize two MUB as long as its phase is 
within [45°,315°]. It is capable of realizing a complete set of MUB if the phase is within 
[111 .5°, 141.7°] or the symmetric interval about 180°. The systematic error in the realization of 
MUB in one wave plate setting is calculated to be twice smaller than that in the conventional 
QWP-HWP setting. TWPs are applied to single-qubit and two-qubit quantum state tomography 
experiments and experimentally show an error reduction by 50% compared with the QWP- 
HWP combination. Other applications of TWP and arbitrary phase plates in the realization of 
any SU(2) and polarization state transformations need to be explored in the future. 

A. Theoretical solutions at intersections 

A.l. Theoretical solutions at the intersection around 8 = 120° 

From the numerical solutions in Fig. [3] some solutions cross at about 8 = 120°. Here solutions 
at this intersection are analytically calculated and the phase at this intersection is indeed 120°. 
We denote the three red solid lines from bottom to top as 0,. i = 1,2,3. The green, blue and 
black solid lines correspond to 90— 0,,9O+ 0, and 180— 0,.i = 1,2,3. Numerical solutions in 
Fig. [3] show that crossing lines at the intersection should have the same coordinates, i.e., 

90° — 0i = 90° + 0i, 02 = 03-90°. (15) 

From Eq. (IT5l >. 6\ = 0°. Substituting 6\ into Eq. (0, one obtains r(0Q = (0,0, l) r . With these 
relations, the conditions that r(0,-),i = 1,2,3 should be orthogonal to each other give two inde¬ 
pendent equations 

cos 5 sin 2 2 02 + cos 2 200 = 0, sin 2 5 sin 2 2 02 = —. (16) 

The solutions of Eq. (IT6l > are 8 = 120°, 02 = ^arccos(—« 27.37° and 03 = 117.37°. Other 
solutions at this intersection represented by green, blue and black can be calculated from their 
symmetries about the red lines. 

A.2. Theoretical solutions at intersections around 8 = 126.3° and 141.7° 

Numerical solutions in Fig.[3]also show some solutions intersect near 8 = 126.3° and 141.7°. 
Here we analytically calculate the solutions and phases at the intersections. The three red lines 
are denoted as 0,, i = 1,2,3 from bottom to top. In Fig.[3]numerical solutions at the intersection 
around 8 = 141.7° give 

e,+0 2 =9O o , 0 3 -90° = 180°-0 3 . (17) 

From Eq. O- 03 = 135° for the intersection around 5 = 141.7°. For the intersection at about 
8 = 126.3°, one has 

0J + 03 =90°, 90°- 0 2 = 0 2 . (18) 

Eq. (fl8l > gives 02 = 45° for the intersection around 8 = 126.3°. With Eqs. ( fiTl i and ( fl8l . con¬ 
ditions of a complete set of MUB at the two intersections give the same equations 

cos 2 8 sin 2 2 0i + cos 8 cos 2 20] — sin 2 8 sin20i = 0, 

(1 — cost)) 2 sin 2 20i cos 2 20i = sin 2 8 sin 2 20i + (cos 8 sin 2 20i + cos 2 20i) 2 . 


( 19 ) 


The two equations in Eq. ( fl9] > are equivalent to 


x(x - 1 )y 2 + x = (1 - x 2 )y, 

(x-l) 2 /= (x-l)V-i 

with x = cos 8 and y = sin 20]. The square of the first equation in Eq. ( l20l > is 
x 1 (x— l) 2 y 4 + 2x 2 (x — l)y 2 +x 2 = (x— l) 2 (x+ 1 ) 2 y 2 ■ 

Replacing (x— 1 ) 2 y 4 in Eq. (OH with the second equation in Eq. Pol l, one obtains 

r 2 

,2 _ x 
y _ 2(1 —x 2 )' 

Substituting Eq. (l22l > into the second equation in Eq. POk we have 

3x 4 + 4x + 2 = 0. 

This four-order equation has two real solutions x = —0.5923 andx = —0.7854, corresponding to 
8 = 126.32° and 141.76°, 0i = 15.66° and 31.90°. FromEq. ®, 0 3 = 74.34° for 8 = 126.32° 
and Eq. ® gives 0 2 = 58.10° for 8 = 141.76°. 

B. Systematic error in the realization of MUB in the QWP-HWP setting 

In the QWP-HWP setting, the realized bases are denoted as r(8 q ,8h,q,h), where 8 cp 5 % are the 
real phases of QWP and HWP in Fig. [I] a). From Eq. © and |[351 , at 8 rl = and 5/, = : r, 

ll r 5,II 2 = sin 2 (4/i-2^)= (r 2 ) 2 , ||r 5 J| 2 = sin 2 2/z, 

ll r c/l| 2 = 4 + 4cos 2 (4/2 — 2q) = 8 — 4(r 2 ) 2 , ||r,|| 2 = 16. 

From Eqs. dlTb and (l24l >. one obtains 

= t ki) 2 = 1. (%) 2 = I 8 - 4(r 2 ) 2 = 20, 

j= 1 7=1 

3 

(e /,) 2 = 48, (e 5 J 2 = ^ sin 2 2 h J . 

7=1 

(eg ; ) 2 depends on the choices of MUB. In the realization of three Pauli operators with /; 1 = 
h 2 = 22.5° and h 3 = 0°, (eg,,) 2 = 1 and 

e 2 = 48(A/i) 2 + 20(A^) 2 + (A5/,) 2 + (A 8 q ) 2 . (25) 

For simplicity, h J is assumed to be uniformly distributed within [0°,360°] in the consideration 
of general MUB. Then (£s,,) 2 i s averaged as |, and 

e 2 = 48(A/7) 2 + 20(A<7) 2 + 1.5(A5/,) 2 + (A8 q ) 2 . 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(26) 



C. Systematic error in qubit state estimation 

From B351I36I I. the systematic error in qubit state estimation with a complete set of MUB is 


tr(P - P ) 2 = \ EII 


( 27 ) 


where ^Jc~' anc * s * s the Bl° c h vector of p. 

In the TWP setting with rotations angles at 0°, 27.37° and 117.37°, from Eqs. (0 and ®, 
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(28) 


Thus, the systematic error tr(p — p) 2 in the estimation of the three states p = p |0) (0| + (1 — 
p)I/ 2 with |0) = | H) , \D) , |/?) is calculated from Eqs. (f27l) and (l28l > to be [8(Af) 2 + |(A S t ) 2 ]p 2 , 
[^(A/) + i(M) 2 ]p 2 and [§(Af) 2 + g(A<5,)]p 2 . 

In the QWP-HWP setting, from (36l . the systematic error tr(p — p) 2 in the estimation of the 
three states can be calculated to be [8(A h) 1 +4(A q) +0.25(A(5/,) 2 ]p 2 , [2(A q) + 0.25 (At)/,) 2 + 
(A 8 c/ ) 2 ]p 2 and [16(A h) 2 +4(A q)\p 2 . 
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